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We present a consistent BPS framework for a generalized Maxwell-Chern-Simons-Higgs model. 
The overall model, including its self-dual potential, depends on three different functions, h (\<p\ , N), 
w (\4>\) and G {\4>\), which are functions of the scalar fields only. The BPS energy is proportional to 
the magnetic flux when w (\4>\) and G (\<j>\) are related to each other by a differential constraint. We 
present an explicit non-standard model and its topologically non-trivial static configurations, which 
are described by the usual radially symmetric profile. Finally, we note that the non-standard results 
behave in a similar way as their standard counterparts, as expected, reinforcing the consistence of 
the overall construction. 



PACS numbers: ll.10.Kk, ll.10.Lm, 11.15. Yc 



I. INTRODUCTION 

Topologically non-trivial structures have been inten- 
sively studied in many areas of physics In particular, 
in the context of classical field theories, these structures 
are described as finite-energy solutions to some non-linear 
models. In this case, such models must be endowed by 
a spontaneous symmetry breaking potential for the mat- 
ter self-interaction, since topological solutions are formed 
during symmetry breaking phase transitions. 

In this context, the most common topological defects 
are kinks Q, vortices [|| and magnetic monopoles Q: 
the first ones are one-dimensional structures described 
by a single real scalar field, while vortices and monopoles 
are two- and three-dimensional configurations arising as 
static solutions of some Abelian- and a non-Abclian- 
Higgs models. 

In particular, vortices are stable configurations de- 
scribed by radially symmetric profile. The simplest 
version of such structures arises as electrically non- 
charged solution of a planar Maxwell-Higgs model en- 
dowed by a fourth-order Higgs potential. However, un- 
der special circumstances, vortices also arise as electri- 
cally charged configurations of some Chern-Simons- and 
Maxwcll-Chcrn-Simons- Higgs (MCS-Higgs) models @. 
In all these cases, such structures can be found as numer- 
ical solutions to a set of first-order differential equations, 
named BogomoFnyi-Prasad-Sommerfield (BPS) ones @. 
In this case, as finite-energy solutions, they have inter- 
esting applications, mainly concerning the superconduc- 
tivity phenomena [3j. 

Moreover, during the last years, beyond the standard 
models cited above, non-usual ones have been intensively 
studied. These theories, generically named k-field models, 
are endowed by non-usual kinetic terms, which change 
the dynamics of the overall model in a non-standard way. 
Here, it is important to point out that the motivation 
regarding such generalization arises in a rather natural 
way, in the context of string theories. 



In fact, k-field theories have been used as effective mod- 
els mainly in Cosmology, as an attempt to explain the 
actual accelerated inflationary phase of the universe Q 
via the so-called k-essence models Q . Furthermore, they 
have been applied in the study of stro ng g ravitational 
waves @ , dark- [l(| and tachyon- matter [lll | and others 
0. 

In such a non-standard scenario, topologically non- 
trivial configurations, named topological k-solutions, can 
exist even in the absence of a symmetry breaking po- 
tential [13j , from which one notes that the existence of 
such structures are quite sensible to the presence of non- 
standard kinetic terms. On the other hand, as an attempt 
to study topological k-solutions via the comparison be- 
tween them and their canonical counterparts, some of 
us have already considered the existence of such solu- 
tions in the context of symmetry breaking k-field mod- 
els [l4j|. Moreover, interesting results concerning these 
models and the corresponding solutions can be found in 
Ref.[H. 

In general, given non-trivial kinetic terms, k-field mod- 
els can be highly non-linear, and the corresponding k- 
solutions can be quite hard to find. In this case, the de- 
velopment of a consistent self-dual theoretical framework 
is quite useful and desirable, since it helps to find topolog- 
ical k-structures as solutions to some non-standard BPS 
equations. Here, as in the usual approach, such equations 
can be obtained via the minimization of the energy func- 
tional related to the non-standard model, the resulting 
self-dual k-solutions being the minimal energy configura- 
tions possible ever. 

In this sense, in a recent work, some of us have pre- 
sented a BPS theoretical framework consistent with a 
generalized self-dual Maxwell-Higgs model endowed by 
non-usual kinetic terms to both gauge and scalar fields 
[Tfjj . Now, we introduce an extension of that work. Here, 
we develop a general first-order approach consistent with 
a non-standard self-dual MCS-Higgs model. The over- 
all model, including its self-dual potential, depends on 
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three different functions, h(\(f>\,N), w (\<p\) and G(|0|), 
which are functions of the scalar sector of the model. In 
this context, we look for topologically non-trivial config- 
urations via the usual radially symmetric static Ansatz, 
and the finite-energy generalized numerical solutions we 
found behave in the same general way as their standard 
counterparts, as expected. 

This letter is outlined as follows: in the next SeclTIl 
we introduce the generalized model and develop the BPS 
framework which allows us to get to its non-standard 
BPS equations. Also, we verify the consistence of the 
overall construction by using it to develop a generaliza- 
tion of the usual MCS-Higgs case, the general model be- 
ing controlled by two real parameters, a and b. Then, in 
Section IIII| we perform the numerical analysis concern- 
ing the non-standard BPS equations previously presented 
by means of the relaxation technique. Also, we depict 
the corresponding self-dual k-solutions for the electro- 
magnetic sector and comment on the main features they 
engender. Finally, in Sec lIVl we present our ending com- 
ments and perspectives. 

From now on, we use the natural units system and a 
plus-minus signature for the planar Minkowski metric. 



II. THE MODEL 

The planar Lagrangian density describing the general- 
ized Maxwcll-Chcrn-Simons-Higgs model is given by 

C d = -m^l F ^^ PAiiFvp (1) 
+w (|0|) \D^\ 2 + hm 2 N) d^Nd»N - U (H , N) , 

where F^ v = d^A v — dyA^ is the usual electromag- 
netic field strength tensor, D M = d^tf) + ieA^cj) is the 
covariant derivative of the Higgs field and & vp is the 
(1 + 2) —dimensional Lcvi-Civita's tensor (with e 012 = 
+1). The additional real scalar field N provides the sta- 
bilization of the self-dual solutions arising in the pres- 
ence of both Maxwell and Chern-Simons terms. Here, 
h (\(f>\ ,N) and w (\4>\) are positive-definite dimensionless 
functions of the scalar fields of the model. The Higgs 
potential U(\(f>\,N) supporting spontaneous symmetry 
breaking is supposed to have the following structure 



u(\<j>\,N) 



fc 2 {N+cmY 

2h(\c/>\,N) 



e 2 A 2 H 2 w (|0|), (2) 



from which one notes that U (\<fi\ ,N) is defined in terms 
of h(\(f>\,N), w (|0|) and G(|0|), with dimG = 1/2. It 
is worthwhile to point out that the specific form of the 
potential in assures the self-duality of the generalized 
model ((T|), i.e., it is a consequence of the BPS construc- 
tion and can be derived from the energy density (jlOl) . 

We introduce the mass scale M of the model, and use 
it to perform the scale transformations: —¥ M~ 1 x fJ ', 



e — > A/ 1 / 2 e and v — s- M x / 2 v, where v is the vacuum 
expectation value of the Higgs field. Then, we get that 
G -> M 1 / 2 G and C d -> M 3 £ , and the modified model 
(fTj) is now described by the dimensionless Lagrange den- 
sity C , which has the same form as Cd- Also, for sim- 
plicity, we choose e = v = k = 1. 

It is well-known that Chern-Simons theories only ex- 
hibit electrically charged static solutions. In fact, the 
static Gauss law related to ([T|) is (j runs over spatial 
indices only) 



d J (hd 1 A°)+2\<f)\ 2 A w = F ] 
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(3) 



from which one observes that, even in the presence of 
non-trivial h(\<j>\ ,N), w (\4>\) and G(|</>|), the absence of 
electrically uncharged static solutions still holds, since 
the temporal gauge (^4° = 0) does not solve ^j. 

Now, since the gauge A = cannot be implemented 
and supposing that h(\<j)\,N), w (\<j>\) and G(|0|) are 
non-trivial functions, the looking for static solutions for 
the second-order Euler-Lagrange equations related to ([1]) 
can be a quite hard task, even in the presence of suitable 
boundary conditions. In this context, the adequate im- 
plementation of the BPS formalism is quite useful, since 
it helps to find finite energy configurations as solutions 
of a given set of first-order differential equations. 

So, from now on, we focus our attention on the de- 
velopment of a BPS framework consistent with (JXJ) and 
@ . Specifically, we look for radially symmetric solutions 
using the standard static Ansatz 

(j> (r, 6)=g (r) e lnB and N (r, 6) = tA° (r) , (4) 



A(r, 



(a (r) -n), 



(5) 



where (r, 9) are polar coordinates and n = ±1, ±2, ±3 ... 
is the winding number (vorticity) of the configuration. 
The fields g(r), a(r) and A (r) must obey the usual 
boundary conditions 

g(0) = 0, o(0) = n, A' (0) = 0, (6) 

3(00) = 1, o(c»)=0, A)M=0, (7) 

where prime denotes the derivative with respect to r. 

Now, we implement the Bogomol'nyi approach on the 
energy functional related to (Q]). Given the generalized 
potential ©, the radially symmetric non-standard en- 
ergy density can be written as 



h fldaV , fdA° 
£ = 2\rTr) +h {-aV 

dg\ 2 , g 2 a 2 ^ 



2h 




2g 2 A 2 Wl 



(8) 



M 1 / 2 (f). N -> M^ 2 N, A p -> M^ 2 A>", k 



Mk, 



where h = h(g,A ), w = w (g) and G = G(g). It is 
clear that, in order to guarantee a positive-definite en- 
ergy, both h and w must be positive-definite and finite. 
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Now, by using the static Gauss law ([3]) together with the 
constraint 



dG 
dg 



2wg, 



(9) 



the energy density ([HJ can be rewritten in the form 

± iA (oG)+ i*( rM °^Y (10) 

r dr r dr \ dr J 

At this point, we see that the resulting total energy is 
minimized by the non-standard first-order equations 



dr r ' 



i^=4 {GtA o ) 

r dr h 



(11) 



(12) 



Equations (|lip and (jT2J) are the BPS equations related 
to the generalized MCS-Higgs model (JXJ) . In the present 
case, the BPS total energy Et ps related to the solutions 
of (jTTJ) and (|12l) can be explicitly evaluated 

E bps = 2n J re hps dr = ±2im, (13) 

where the minimal energy density eb ps is defined as 

, 1 d , ^ 1 d ( , *a dA °\ 
e bps = ±-—(aG) + -—(rhA —). 14 

r dr r dr \ dr J 



r dr 



Also, to compute the BPS energy (|T3)) . we have used the 
boundary conditions ([6]) and ([7|) together with 

G(r = 0) = -1, G(r = oo) = 0. (15) 

Finally, we assume that /i (\cf>\ ,N) and u; (|</>|) behave as 

h(r = 0) = H , h(r = oo) = Hoo, (16) 



w (r = 0) = Wq, w (r = oo) = Wo, 



(17) 



where i?0; Wo an( i ^oo are real non- negative constants, 
and Woo is a real positive constant. 

From p3|) . one notes that the total energy Eb ps is 
quantized according the winding number n. Further- 
more, it can be related to the magnetic flux $b in the 
standard way. Also, we point out that the first-order 
framework developed here is implemented for any func- 
tion h{\4>\ , N) finite and positive-definite, i.e., given some 
function h(\<j)\,N), the set formed by (fTTj) and (|T2"j) gives 
the BPS configurations of the generalized MCS-Higgs 
model ([TJ. On the other hand, in the absence of ([2]) 
or ([5]), the energy functional (JSJ cannot be written as 
(|10p. and the development of a consistent non-standard 
first-order formalism cannot be performed. 



The prescription for the development of generalized 
self-dual Maxwell-Chern-Simons-Higgs models is as fol- 
lows: given any function w (\<j>\), positive and finite, the 
corresponding G(|</>|) is obtained by means of Eq. ©. 
Then, given any function h(\<p\,N), also positive and 
finite, the resulting self-dual potential U(\<fi\ ,N) is de- 
termined via ©. In this scenario, the BPS states of the 
generalized model <m are the solutions of the equations 
©, (JTT) and (HI]). The total energy is given by and 
the energy density by (fTl)) . 

In the standard MCS-Higgs case, BPS solutions with 
non-trivial topology only exist in the asymmetric vacuum 
of the potential 



U s (\<P\,N) = ^(\0\ 2 +N-iy + N 2 



(18) 



which is defined by N = and \cj>\ = 1. So, for sim- 
plicity, we assume that also the generalized potential $2$ 
achieves its asymmetric vacuum when N = and \<f>\ = 1. 
As a consequence, topologically non-trivial self-dual so- 
lutions of ([T]) exist in the same asymmetric phase as their 
standard counterparts. 

The standard MCS-Higgs model is trivially recov- 
ered starting from our generalized framework by setting 
w(|^|) = h(\4>\,N) = 1. So, instead of recovering the 
usual model, we introduce a generalization of such the- 
ory. The generalized model is defined by 



w (H) = &(|0| 2 -i) 



6-1 



(19) 



where b is a positive odd- number. The corresponding 
G (|0|) is 

G(|0|) = (V| 2 -l)\ (20) 
For simplicity, we choose h(\cj>\ ,N) as follows 



h(N) = aN 2 + 1, 



(21) 



where a is a real non-negative number. Note that a = 
and 6=1 leads us back to the standard case. The 
resulting generalized self-dual potential U (\4>\ ,N) ^ is 



U(\</>\,N) = 



N ■ 



(H 2 -i)' 



2(cdV 2 + i) 



+WV 2 |0| 2 (|0| 2 -l) b \ (22) 



In the present case, <j3j> , (fTTj) and ([T2|) can be written, 
respectively, as 



i ,2 ,n fd 2 A° ldA°\ n A0 fdA°\ 



dr 2 r dr 
2bg 2 (g 2 -!)"- 



i\b-i a0 _ Ida 
r dr 



(23) 
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dg = ag 
dr r ' 



Ida _ ±(.9 2 -l) - A 
r dr 



1 



(24) 



(25) 



which must be solved according the finite energy bound- 
ary conditions ([6]) and ([7]). 

In the next Section, we solve ([23]) . (f24"]) and ([25]) via © 
and ((?]) for different values of a and 6. Then, we use the 
numerical solutions we found for g (r), a (r) and A (r) 
to depict the corresponding profiles for the electric field 



the magnetic field 



dA° 



Ida 

B (r) = — — , 
r dr 



(26) 



(27) 



and the minimal energy density (|14j) . Also, we comment 
on the main features the non-standard numerical solu- 
tions engender. 



E(r) 




FIG. 1: Solutions to E (r) for a — and 6=1 (dash-dotted 
black line), a = and b = 3 (dotted blue line), a — 5 and 
6=1 (dashed red line) and a = 5 and 6 = 3 (solid green line). 



III. NUMERICAL SOLUTIONS 

In this Section, we focus our attention on the non- 
standard numerical solutions themselves. The equations 
to be studied arc ([23]). (|24| and ([25} . and the fields g (r), 
a (r) and A (r) must behave according the finite energy 
boundary conditions ([6]) and ([7]). 

Here, for simplicity, we choose n = 1. Then, we solve 
the first-order system numerically by means of the relax- 
ation technique, for different values of the real parameters 
a and b. The profiles of the physical fields and energy 
density are depicted in the figures below. The system 
was solved for a = and 6=1 (usual case, dash-dotted 
black line), a = and 6 = 3 (dotted blue line), a = 5 
and 6=1 (dashed red line) and a = 5 and 6 = 3 (solid 
green line). 

In Figure 1, we depict the numerical solutions for the 
electric field ([26]) , and we see that the non-standard ones 
behave, in general, as their usual counterpart: starting 
from (zero), the solutions reach their maximum val- 
ues at some finite distance R from the origin, and vanish 
as r goes to infinity. However, we point out that differ- 
ent solutions exhibit not only different amplitudes, but 
also slightly different characteristic lengths: in general, 
increments on a and/or 6 lead to decrements on ampli- 
tudes and/or characteristic length. Even in this case, 
such variations are expected to occur in the context of a 
non-standard theory as ([1]); see, for instance, fl4| . 

In Figure 2, we present the profile of the magnetic field 
(|27|). The same manner as for the electric field, the gener- 
alized solutions behave as the standard one: the magnetic 
fields arc lumps centered at the origin, and they decrease 
monotonically as r goes to infinity. In this case, for a fixed 




r 



FIG. 2: Solutions to B (r). Conventions as in FIG. 1. 

a and an increasing 6, the amplitudes experience a small 
increase and the characteristic length decreases. On the 
other hand, for a fixed 6 and an increasing a, the am- 
plitude will decrease significantly, but the characteristic 
lengths remain approximately fixed by a compensatory 
effect relating a and 6. 

Numerical solutions for the minimal energy density 
£b P s (ff4]) are depicted in Figure 3, and we see that also 
such solutions are lumps centered at r = which decrease 
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FIG. 3: Solutions to e 6ps . Conventions as in FIG. 1. 

monotonically to r — > 0. The lumps profile are similar 
to those the magnetic field. However, in this case, for a 
fixed a and an increasing b, the amplitudes experience 
a large increase and, the characteristic lengths decreased 
slightly. On the other hand, for a fixed b and an increas- 
ing a, the amplitude will decrease significantly but the 
characteristic lengths remain almost fixed by a compen- 
satory effect relating a and b. 




r 



FIG. 4: Solutions to B (r) for n = 5. Conventions as in FIG. 
1. 

To end this Section, we have depicted the numerical 



profiles for B (r) and £b ps for n = 5 (see figures 4 and 
5, respectively). It is well known that in the standard 
model, for a large winding number, the profiles of both 
the magnetic field and the energy density go from a lump 
(n = 1) to a ring (n >> 1). From figures 4 and 5, one 
notes that also in this limit (increasing vorticity), the 
generalized solutions mimic the usual ones, as expected. 
For large n with fixed a and a increasing b, the profiles 
experience a slight change. However for b fixed and an 
increment in a, the profile's amplitudes experience a con- 
siderable decrease. 



IV. ENDING COMMENTS 

In the present letter, we have performed the develop- 
ment of a first-order theoretical framework in the context 
of a generalized MCS-Higgs model given by ([T]). In the 
present case, the non-standard model, including its self- 
dual potential ([2]), is given in terms of three different 
functions, h (\(j>\ , N), w (\(j>\) and G (\</>\), which are func- 
tions of the scalar fields only. Here, in order to avoid 
problems with the energy of the model, h(\<j>\ ,N) and 
w(|(/)|) must be positive-definite. Then, given the gen- 
eral structure for the self-dual potential @, the consis- 
tence of the non-standard first-order approach only holds 
when u> (|</>|) and G(|<^|) are related by the differential 
constraint ©. On the other hand, one notes that there 
is no additional constraint to be imposed on h (|</>| , N). 
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FIG. 5: Solutions to Eb ps for n = 5. Conventions as in FIG. 
1. 

After perform the construction of the generalized 
MCS-Higss model, we illustrate such realization provid- 
ing an explicit non-standard model specified by (jlU)) , ([20)1 
and pTj) . In this case, the non-standard model is con- 
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trolled by two real parameters, a and b. Immediately, we 
have investigated radially symmetric self-dual configura- 
tions given by the usual static Ansatz ((4]) and (0, and 
considering that the fields behave according the usual fi- 
nite energy boundary conditions (O and 0. 

We have integrated the generalized BPS equations by 
means of the relaxation technique, for different values of 
a and b. The numerical profiles we found are depicted 
in figures 1, 2, 3, 4 and 5. In general, we have seen 
that the non-standard solutions mimic their usual coun- 
terparts. Also, as expected, we have noted variations on 
the amplitudes and on the characteristic lengths of the 
modified solutions. Furthermore, we have identified the 
rules controlling such variations. 

This work is an extension of a recent investigation per- 
formed by some of us [l6|. In this sense, a very interest- 



ing issue to be considered concerns the development of a 
non-standard self-dual MCS-Higgs model which exhibits 
the very same numerical solutions engendered by the 
standard theory. In this case, such investigation would 
be a generalization of the one presented in [l7j , regard- 
ing a non-standard Maxwcll-Higgs model. Another issue 
concerns the supersymmetric extension of the self-dual 
model (H)-©. These issues are in advance, and we hope 
report interesting results in a near future. 
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